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Abstract. The results from hydrodynamical TREESPH simulations of 
galaxy clusters are used to investigate the dependence of the final cluster 
X-ray properties upon the numerical resolution and the assumed star 
formation models for the cooled gas. When cold gas particles are allowed 
to convert into stars the final gas profiles show a well defined core radius 
and the temperature profiles are nearly flat. Star formation methods 
based on a local density threshold, as in Navarro and White (1993), are 
shown to give stable results. Final X-ray luminosities are found to be 
numerically stable, with uncertainties of a factor ~ 2. 



1. Introduction 

Hydrodynamical simulations have been widely used to predict for different the- 
oretical frameworks the time evolution of the temperature distribution of X-ray 
clusters. With increasing availability of computational power numerical hydro- 
dynamical simulations have been attempted to model the effects of radiative 
cooling on the gas in the formation and evolution of cluster galaxies (Katz & 
White 1993; Anninos & Norman 1996; Yoshikawa, Jing & Suto 2000; Pearce et 
al. 2000; Lewis et al. 2000). In this contribution I will show the preliminary 
results that have been obtained from a series of hydrodynamical SPH simula- 
tions of galaxy clusters. The simulations have different numerical resolution 
and include the effects on the gas component of radiative losses, star formation 
and energy feedback from SN. Final profiles are compared in order to assess the 
effects of numerical resolution, or different star formation prescriptions, on the 
cluster X-ray variables. 



2. Simulations 

In a previous paper (Valdarnini et al. 1999, hereafter VGB) a large set of hydro- 
dynamical simulations was used to study the global X-ray cluster morphology 
and its evolution. The simulations were run using a TREESPH code with no gas 
cooling or heating. I refer to VGB for a detailed description of the simulations. 
In Valdarnini (2001) I have used the same cluster sample to study the effects 
of including in the simulations additional physics such as gas cooling and star 
formation. I will report here the results for the cluster with label ACDMOO in 
VGB. This is the most massive cluster (M^ ~ 1.5 • W^^Mq) extracted from a 
cosmological ACDM N-body simulation with size L = 200h~^ Mpc, matter den- 
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sity = 0.3 and Hubble constant Hq = 70Km sec Mpc . In order to check 
the effects of radiative coohng for this cluster a set of TREESPH simulations 
was performed, with initial conditions provided by the cosmological simulation. 
The effects of radiative cooling are modelled in these simulations by adding to 
the SPH thermal energy equation an energy-sink term represented by the radia- 
tive cooling function. The simulation parameters are: number of gas particles, 
Ng 22, 500; gas particle softening parameter £g ~ 30Kpc; gas particle mass, 
nig ~ 10^*^ Mq; number of dark particles, Nd ~ 45,000 (including those in an 
external shell surrounding the cluster, cf. VGB). 

Allowing the gas to cool radiatively will produce dense clumps of gas at low 
temperatures 10^ °K). In these regions the gas will be thermally unstable 
and will likely meet the physical conditions to form stars. In TREESPH sim- 
ulations star formation (SF) processes have been implemented using different 
approaches . According to Katz, Weinberg & Hernquist (1996, hereafter KWH) 
a gas particle is in a star forming region if the flow is convergent and the local 
sound crossing time is larger than the dynamical time (i.e. the fluid is Jeans un- 
stable). In a simplified version Navarro & White (1993 , hereafter NW) assume 
as a sufficient condition that a gas particle must be in a convergent flow and 
its density exceeds a threshold {pg > pg^c = 7 ■ 10~^^ ^rcm"^). If a gas particle 
meets these criteria then it is selected as an eligible particle to form stars. The 
local star formation rate (SFR) obeys the equation 

dpg/dt = -C^Pg/Tg = -(Ip^/dt , (l) 

where pg is the gas density, Pi, is the star density , c^, is a characteristic dimen- 
sionless efficiency parameter, Tg is the local collapse time and it is the maximum 
of the local cooling time Tc and the dynamical time r^. At each step the prob- 
ability that a gas particle will form stars in a time step At is compared with a 
uniform random number. If the test is successful then a mass fraction of the 
gas mass is converted into a new coUisionless particle. This star particle has the 
position, velocity and gravitational softening of the original gas particle. Typical 
assumed values are = 1/3 and = 0.1 (KWH). For the local SFR, NW adopt 
Eq. 1 with Ci, = 1, Tg = Td and = 1/2 when a gas particle can convert part of 
its mass into a star particle. The numerical tests have been performed following 
the NW prescriptions for selecting gas particles which can form stars. Once a 
star particle is created it can release energy into the surrounding gas through 
supernova (SN) explosions. All the stars with mass above 8Mq end as a SN, 
leaving a IAMq remnant. The SN energy (~ lO^^er^r) is released gradually into 
the gas according to the lifetime of stars of different masses. 

3. Results 

Three simulations were performed with the same numerical parameters, but 
using different SF methods or parameters. This in order to test the final de- 
pendence of the cluster X-ray variables on the assumed star formation model 
for the cooled gas. In simulation I cold gas particles have been converted into 
stars according to the NW method. In the other two runs (II and HI) the KWH 
prescription is adopted for converting gas particles into stars, but with different 
values of the star formation efficiency parameter c*: 0.1 and 1.. As is shown in 
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Figure 1. Plots showing several cluster properties in simulation 

runs with different SF prescriptions. The continuous line refers to the 
NW method, the others to KWH with different (c in the bottom left 
panel). Top left: ratio of the star mass within the radius r over the gas 
mass within r. Top right: final radial density behavior for the gas com- 
ponent. The simulation results are compared with the corresponding 
high-resolution runs (see text). To facilitate a comparison the radial 
profiles of the high resolution results have been shifted downward by 
lO'^, /c = 2 for densities and A; = 1 for temperatures. Bottom right: Ra- 
dial temperature profiles. Bottom left : X-ray luminosity versus time, 
the thick lines correspond to the high resolution simulations. 
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Fig. 1 the most important differences are found between the KWH simulations 
with different c*. The differences are dramatic in the final X-ray luminosities, 
which differ by a factor ~ 40. The source of this discrepancy relies in the dif- 
ferent gas density profiles, which have substantial differences in the cluster core 
regions for r ^ lOOKpc. These differences are localized at the cluster center, 
beyond r ~ WOKpc all the profiles converge, as it is shown in the plots of Fig. 
1. The temperature profiles have a peak value of ~ 10^ °K at r lOOKpc and 
thereafter decline outward by a factor ~ 2 out to r^. Below ~ lOOKpc the 
profiles instead show large differences. Compared to the NW run the simulation 
with c^, = 0.1 has gas temperatures which decline inwards by a factor ~ 10 from 
~ lOOKpc down to r ~ lOKpc. These radial decays follow because of the less 
efficient conversion of the cooled gas into stars compared to the NW run. There 
is a remarkable agreement for the ratio of cluster mass locked into stars to the 
gas mass, which is ~ 10% at for all the runs considered. In order to assess the 
effects of numerical resolution upon final results simulations I, II and III have 
been run again but with a number of particles increased by a factor ~ 3. These 
simulations will be referred as IH, IIH and IIIH, respectively. The simulation 
results are shown in Fig. 1. For simulations IH there are not appreciable differ- 
ences in the radial profiles. The profiles of simulation IIH are instead different 
from those of nin II at r bOKpc. The strong drop in T{r) has been removed 
and the gas density profile is much closer to the NW one. Simulation IIIH yields 
final profiles very similar to the ones of the parent simulation. The bottom left 
panel of Fig. 1 shows that high resolution runs have final X-ray luminosities 
which can differ within a factor ~ 2 from the parent simulations. 

To summarize, the above results demonstrate that simulations I and HI have 
an adequate numerical resolution to reliably predict X-ray cluster properties, 
such as the X-ray luminosity. For simulation II ( KWH with = 0.1 ) there are 
large differences at the cluster core between the final profiles when the numerical 
resolution is increased. 
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